NACA TN 2621 


\ 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 2621 


DEFLECTION AND STRESS ANALYSIS OF THIN SOLID WINGS 

OF ARBITRARY PLAN FORM WITH PARTICULAR 

REFERENCE TO DELTA WINGS 

By Manuel Stein, J. Edward Anderson, 
and John M. Hedgepeth 

Langley Aeronautical Laboratory 
Langley Field, Va. 

CASE FILE 
COPY 



Washington 
February 1952 


« 


* 


IT 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2621 


DEFLECTION AND STRESS ANALYSIS OF THIN SOLID WINGS 

OF ARBITRARY PLAN FORM WITH PARTICULAR 

REFERENCE TO DELTA WINGS 

By Manuel Stein, J. Edward Anderson, 
and John M.. Hedgepeth 


'SUMMARY 


The structural analysis of arbitrary solid cantilever wings by 
small-deflection thin-plate theory is reduced to the solution of linear 
ordinary differential equations by the assumption that the chordwise 
deflections at any spanwise station may be expressed in the form of a 
power series in which the coefficients are functions of the spanwise 
coordinate. If the series is limited to the first two and three terms 
(that is, if linear and parabolic chordwise deflections, respectively, 
are assumed) , the differential equations for the coefficients are solved 
exactly for uniformly loaded solid delta wings of constant thickness and 
of diamond chordwise cross section with constant thickness ratio. For 
cases for which exact solutions to the differential equations c ann ot be 
obtained, a numerical procedure is derived. Experimental deflection 
and stress data for constant-thickness delta-plate specimens of V?° 
and 60 ° sweep are presented and are found to compare favorably with the 
present theory. • 


INTRODUCTION 


One of the present trends in the development of high-speed air- 
planes and missiles is toward the use of thin low-aspect-ratio wings. 

The structural analysis of these wings often cannot be based on beam 
theory since the structural deformations may vary considerably from 
those of a beam and, indeed, may more closely approach those of a plate. 
In cases where the wing construction is solid or nearly solid the use of 
plate theory in the analysis is particularly valid, and it is this type 
of wing which is considered in the present paper. 

Exact solutions to the partial-differential equation of plate 
theory are not readily obtained, especially for plates of arbitrary 
shape and loading; however, a number of approximate solutions to 
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specific problems on cantilever plates have appeared in the literature 
(see, for example, references 1 to 7) • Of the approaches used in these 
references, only the one in references 6 and 7 is readily applicable to 
plates of arbitrary plan form, thickness distribution, and load distri- 
bution; thus it is the most useful one for the analysis of actual wings. 


In reference 6 the cantilever-plate problem is simplified by the 
assumption that the deformations of the plate in the chordwise direction 
(parallel to the root) are linear. By minimizing the potential energy 
of the plate, the. partial-differential equation of plate theory is 
replaced by two simultaneous ordinary differential equations for the 
spanwise variations of the bending deflection and twist. In reference 7 
the same ordinary differential equations are obtained in a different 
manner. Refinement of the analysis by inclusion of the effect of 
parabolic, cubic, or higher-order chordwise camber terms is indicated 
in reference 6, and as the order of refinement is. increased a corre- 
sponding increase in the number of ordinary differential equations is 
obtained. 

In the present paper, which is an extension of reference 6, a general 
set of ordinary differential equations is presented which may be used to 
obtain. any desired degree of approximation to the deflection of the plate. 
These equations are solved exactly for several cases of delta plates under- 
uniform load first by considering linear chordwise deformation only and 
second by including the effect of parabolic chordwise camber. Comparisons 
are drawn between the stresses and deflections computed from the equations 
of each approximation and also with some experimental results. 

The differential equations presented contain coefficients that 
depend on -the plan form and stiffness distribution of the plate and on 
the loading. In the present paper, the plates considered in detail have 
coefficients such that the differential equations can be solved exactly; 
however, in cases for which exact solutions cannot be obtained a numeri- 
cal procedure must be used. One such procedure is derived and its 
accuracy is demonstrated. 


SYMBOLS 


l length of plate measured perpendicular to root 

c root chord of plate 

lateral load per unit area, positive in z-direction 


P 
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D 

E 


local thickness of plate 
average thickness of plate 

local flexural stiffness [■ — r 

H 1 - “ ) 

flexural stiffness based on average thickness 

modulus of elasticity of material 
Poisson* s ratio 




w deflection of plate, positive in z-direction 

x>y,z coordinates defined in figure 1 

^n function of x, coefficient in power series for deflection 

w 

n=0 

c l(x), c 2(x) functions defining plan form (see fig. l) 


RESULTS 


variable obtained by transformation xi. = 


a x>°y 


'xy 


normal stresses 
shear stress 

maximum principal stress 
aspect-ratio parameter I 


'iJSTT 

c][2 


4 ) 


The derivation of the general set of ordinary differential equa- 
tions is given in appendix A. The general procedure outlined in- refer- 
ence 6 is followed; that is, the deflection of the plate w is 
expanded into a power series in y the chordwise coordinate, with 
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coefficients which are functions of x the spanwise coordinate (see 
fig. 1) 

w = qp Q (x) + cp x (x)y + cp 2 (x) y 2 + . . . + cp^x)^ (l) 

Equation (l) is substituted into the expression for the potential energy 
of the plate-load combination which is in turn minimized by the calculus 
of variations with respect to each of the coefficients cp^. The results 
of the variational procedure appear as N + 1 simultaneous differential 
equations with the coefficients qpn as unknowns. 

• By taking a sufficient number of terms in the expansion of w, the 
resulting differential equations can be used to obtain any desired 
degree of accuracy in the solution for the deflections of any given 
cantilever plate subjected to an arbitrary lateral load. Of most 
interest, perhaps, are the particular cases for N = 1 and N = 2, 
which are obtained from the general set of equations and are simplified 
in appendix A. The case for N = 1 (also derived in references 6 and J) 
includes linear chordwise deflections, and the case for N = 2 takes - 
into account parabolic chordwise curvature. Although for most practical 
problems the solution by the parabolic theory should be adequate, cases 
might exist in which cubic, quartic, or even higher-order chordwise 
terms should be included, depending on the convergence of the series for 
the configuration considered. 

The particular equations for N = 1 and N =2 are used to 
determine the deflections and stresses of the following cantilever 
plates subjected to uniform lateral load: 

(1) A 45° delta plate of uniform thickness 

(2) A 60° delta plate of uniform thickness 

(3) A 45° delta plate of diamond chordwise cross section with 

constant thickness ratio 

Fortunately, for these configurations, the solution can be carried 
out exactly by both the linear and parabolic theories, and the details 
of these exact solutions are included in appendix B. In general, how- - 
ever, exact solutions cannot be obtained and some numerical method. must 
be used. One such method, based on replacing derivatives by their 
first-order-approximation difference forms, is derived in appendix C. 

A summary of the results for the three particular problems is shown in 
figures 2 to 11. Deflections obtained by the linear theory and the para— 
'bolic theory for the three configurations are compared in figures 2, 3> and 4 
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The stress results for the 45 and 6o° uniform- thickness delta 
plates indicate that both the linear and the parabolic theories are 
adequate and that the parabolic theory is better than the linear theory 
only near the root. It should be noted that, although the maximum 
principal stress over a large part of the 45° plate is plotted in fig- 
ure 5> only the stresses normal to the root along the line ^ = O.OO 87 

of the 60 ° plate are plotted in figure 6 since only these stresses are 
given in reference 8 . (The maximum principal stress and the stress 
normal to the root are theoretically equal at the root since the root 
shear stress is zero.) 

Experimental data are lacking for the diamond-cross-section delta 
plate and, therefore, only theoretical stresses are shown in figure 7 . 
As in the case of deflections, the results obtained from the linear 
theory and those obtained from the parabolic theory are almost coinci- 
dent, the difference being greatest near the root. Figure 7 has also 
been plotted so that the results are directly comparable with those for 
the 45 uniform- thickness plate in figure 5 on an equal-weight basis. 

As can be expected, the diamond-cross-section, constant-thickness-ratio 
plate is a better design structurally; the stresses in the diamond- 
cross- section plate are everywhere smaller and are almost constant in 
the spanwise direction. 

The theoretical results in figures 2 to 7 have been obtained from 
exact solutions of the differential equations of the linear and 
parabolic theories. In order to test the reliability of the numerical 
method derived in appendix C, the differential equations were also 
solved numerically. The results shown in figures 8 and 9 indicate that 
the agreement is good between the numerical solution in which five equal 
intervals were used and the exact solution of the differential equations 
for the case of the 45° uniform- thickness plate. The same good agree- 
ment can be expected in other cases where the thickness and load dis- 
tributions are not too erratic and where the plate stiffness does not 
go to zero at the tip; that is, when no signularities appear at the tip. 

Since the efficacy of the numerical method depends on how well 
parabolic arcs fit the various functions between stations, serious 
error can result from blind application. An example of the seriousness 
of these errors and of the manner in which they can be remedied is 
shown in figures 10 and 11 . In these figures a comparison is made 
between exact and numerical results obtained on the 45° diamond- cross- 
section, constant- thickness-ratio plate. As can be expected, the five- 
station numerical solution fails to follow the exact solution in the 
neighborhood of the singularity at the tip. Since the region of trouble 
is localized at the tip, a reasonable remedy would be to decrease the 
spacing of the station points near the tip. This decrease in spacing 
may be accomplished either by using a greater number of equally spaced 
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Stresses obtained by the linear theory and the parabolic theory for the 
three configurations are compared in figures 5, 6 , and 7 . Where avail- 
able, experimental deflections and stresses are also shown in these 
figures. The details of the procedure used to obtain the experimental 
deflections of the 45° and 60 ° uniform- thickness plates and the experi- 
mental stresses in the 45° uniform-thickness plate are contained in 
appendix Dj whereas the experimental root stresses for the 6o° uniform- 
thickness plate were obtained from reference 8 . Figures 8 to 11 present 
the comparison between deflections and stresses computed from the exact 
solutions of the differential equations and those computed from the 
numerical solution of the same equations. 


DISCUSSION 


The results shown in figures 2 and 3 indicate that, with regard to 
deflections, either the linear theory or the parabolic theory is 
adequate for the case of a constant- thickness delta plate subjected to 
a uniform load, the comparison being somewhat better for the 60 ° plate 
than for the 4-5 plate. If accurate slopes in the chordwise direction 
(angle of attack) -are desired, however, the parabolic theory must be 
used because the error in the angle of attack as computed 'by the linear 
theory is as much as 30 percent (see figs. 2 and 3) • The appreciable 
anticlastic curvature, evidenced by the experimental results of figures 2 
and 3 , may be important aerodynamically and is, of course, not taken 
into account by the linear theory. 

The apparent convergence of the aforementioned series in the case 
of the diamond-cross-section plate (see fig. 4) implies that the linear 
theory is adequate for this case. The lack of chordwise curvature in 
the result obtained by the parabolic theory is attributable to the fact 
that the natural tendency of the plate to have anticlastic curvature is 
canceled by the opposite tendency of the thin edges to bend down under 
the load. Unfortunately, no experimental results are available for this 
configuration. 

In figure 4 the plate stiffness D in the nondimen sional param- 
eter wD/pZ is the local value of D at a point where the thickness is 
equal to the average thickness of the plate as a whole. Thus the 
results of figure 4 are comparable with the results of figure 2 on an 
equal-weight basis. It can be seen that the deflections of the diamond- 
cross-section, constant- thickness-ratio plate are everywhere less than 
those of the uniform- thickness plate although they increase rapidly near 
the tip. This curling-up or singularity in slope at the tip is a result 
of using a small-deflection theory and probably would not be so marked 
in an actual case. 
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stations or by using unequally spaced stations crowded near the tip. 
The increase in accuracy obtained by increasing the number of equally 
spaced station points to ten is shown in figures 10 and 11. ' 


CONCLUDING REMARKS 


The general method presented herein for finding deflections and 
stresses of solid or nearly solid wings is, in principle, capable of 
yie ing arbitrarily accurate results for any configuration. It is 
seen that, for the examples considered, only the first two or three 
terms in the series expansion need be considered to obtain adequate 

f In f^ ltlon ' f or most Practical plate-like wings with clamped 
roots the first two or three terms will probably be adequate, although 
problems may exist wherein more terms are needed. 

The numerical procedure, derived for application in cases where 
exact solutions cannot be obtained, gives good agreement when' compared 
with exact solutions if enough stations are taken along the span. The 
necessary number of stations is dependent on the type of thickness and 
loading distribution considered, five equally spaced stations being 
enough for . the uniform-thickness delta wing subjected to uniform loading, 

^ in§ - neCeS K ary /r the ^ amon< ^ _cr oss- section, constant-thickness- 
ratio delta wing subjected to uniform loading. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., November 30 , 1951 
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APPENDIX A 

DERIVATION OF DIFFERENTIAL EQUATIONS 


The structure considered herein is a thin, elastic, isotropic, 
cantilever plate of arbitrary plan form and slowly varying thickness 
subjected to distributed lateral load (see fig. l) . By assuming that 
the deflection of the plate can be represented by a power series in the 
chordwise coordinate and by applying the minimum-potential-energy 
principle, a set of ordinary differential equations in the spanwise 
coordinate is obtained from which the coefficients of the power series 
may be determined. First the general set of equations is derived; then 
the particular equations for the cases of linear chordwise deflections 
and parabolic chordwise deflections are deduced from the general set 
and simplified. 

General equations .- The potential energy of the system under 
consideration is 


Potential energy = 


rl p c 2 (x) 
0 ' ^c x (x) 


D(x,y) 


+ (¥*) + 2u & tl + 


2(1 - h) 


/ 

y) 


in which 


D(x,y) = 


E 


12 


,8x 

- p(x,y)w 

[^(x,yj] 3 

(l - p 2 ) 


W7 

dy dx 


8x^ 8y^ 


(Al) 


and p(x,y) is the distributed lateral load. 

The assumption is made that the deflection w can be represented 
by the power series 


w 



n=0 


<P n U)y n 


(A2) 


Substitution of this expression for 


w into equation (Al) gives 
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Potential energy = / dx ^ 

J 0 


1 r " " 

2 A__ a m+n+l < Pm ^n + 
c m=0 n^(3 L 


mn(m - l)(n - Da^^qy^ + 2pn(n - + 


in which 


a r = 


rc 2 (x) 
C]_(x) 

>C 2 (x) 


, n J- 

2(1 - n)mna ^ .,® ® - > p © 

m+n-lMn T n Z_ *n+l T n 
-> n=0 


I* 1 

C(x,y)y dy (r = 1, 2, . . . 2N + l) 


P r = / P(x,y)y r_1 dy (r = 1, 2, . . . N + l) 

J c 1 (x) 


(A3) 


(A4) 




and the primes denote differentiation with respect to x. 

Minimization of the potential energy by means of the calculus of 
variations gives 


6(Potential energy) = 0 


dx^ 


m=0 n=0 L 


+ <P n " 8 0 + ' 


mn(m - l)(n - l)a m+n _ 3 (q^ l &<P n + q> n 6q^) + 
2pn(n - + cp n &<') + 

2(1 - M)rona m + n _ 1 ^9 m , 5(p n ' + <P n 'sO] - 
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Integrating by parts and collecting terms results in 


pi N 

0 = / dx ^ &cp n - 


n=0 


XI K+n+lV) + ^ m ( m - 1) K+n-l^m)" 

m=0 L 

2(1 - 4) m n^a m+n _ 1 q5 m ' )' + pn(n - + 

r 

mn(m - l)(n - l)a m+n _3q\ n - P n+] > + 4X1 &9n' 

J n=0 


U-- 

a m+n+l < Pm + 


um(m - D^^cpjl — fl 8cp n j^n+n+i^")' + 

J o & 7 

pm(m - l)(a m+n . A )’ - 2(1 - u)mna m+n . A '] j 


(A5) 


• Everywhere in the region of the plate, except at the boundary x = 
the variation of w is arbitrary. At x = 0 the cantilever boundary 
conditions 


*-s-° 

yield 

<P n (°) = 9 n '(0) =0 (n * 0, 1, ... N). (A 6) 

and therefore the variation in these quantities must also he zero 
at x = 0. 

Equation (A5) is then satisfied if , in addition to equation (A 6) , 


N 

XI (am+n+lO + - 1 )( a m+n-l c Pm) “ 2(1 - ) + 

pn(n - l) a m+n _ + nm(m - l)(n - l)a m+n _ 3 <P^j = ?r 


' n+1 


(n = 0, 1, . . . N) (A7) 


ki+n+A" + Mm(m - • = 

m=0L J x=Z 


= 0 (n = 0, 1, . . . N) (A8) 
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and 


XI iK+n+W) + Win - l) ( a m+ n - 1 c Pm )' " 2 (l - ^am+n-lVl = 0 

m=0 L Jx=Z 

(n = 0, 1, . . . N) (A9) 

Equations (A7) form a set of N + 1 simultaneous ordinary differ- 
ential equations for the functions cp n (x) . The functions <p n are com- 
pletely determined hy these differential equations and the boundary 
conditions (A 6), (a 8), and (A9) • 

Particular case of linear chordwise deflections .- If N = 1, the 
deflection function becomes 


w = <Pq + y9 x , (A10) 

a linear function in the chordwise direction, where cpQ is the bending 
deflection and is the twist. Equations (A7) become 

( a ]V) + (^l") = Pi (All) 

( a 2V) + ( a 3 9 l") - 2(1 - n^a^' )' = p 2 (A12) 

The root boundary conditions, given by equation (A 6), become 

<P 0 (o) = V (0) = Cp i(°) = = 0 (A13) 

The tip boundary conditions, given by equations (A8) and (A9) , become 


(Vo" + Vi') x=! 

^ a 2^o" + a 3^l"^ ^ 


= 0 


= 0 


(Al4) 


(A15) 
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j( a 2 < Po") ' + - 2(1 - n)a 1 (p 1 |j _^ = 0 (A17) 

Equations (All) to (A17) are the differential equations and corresponding 
boundary conditions presented in reference 6 (if only distributed load 
is considered) where the symbols W and 6 are used instead of qp 0 
and cp}_, respectively. 

If equation (All) is integrated twice and the boundary condi- 
tions (Al4) and (Al6) are used, 



(A18) 


Substitution of cp 0 into equations (A12) , (A15) , and (A17) gives 


’Vl") - 2(1 - ^(a^') =P 2 - 

/ a p r l r l 2 \" 

— / / P 1 dx ] 
\ a l J x J x 1 ) 

(A19) 

/ »\ 

Vl _ = 

0 

(A20) 

\ /x=Z 

[Vi”)' - 2(1 - 

■] 

(A21) 


x=2 


in which 



If equation (A19) is integrated once and the boundary condition (A2l) 
is used. 



- 2(1 - nja^' 



(A22) 


The differential equation (A22) is a second-order differential equation 
in qp-^ . The twist (p-^ and then the bending deflection (p 0 are 

obtained by solving equations (A22) and (Al8) , respectively, by 
applying the boundary conditions (A13) and (A20) . 
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Particular case of parabolic chordwise deflections .- The effect of 
parabolic chordwise camber may be included by letting N = 2 in the 
general power series (equation (A2)). If N = 2 the deflection 
function becomes 


w = 9 0 + y ( P 1 + y qp 2 

Here represents the spanwise distribution of parabolic chordwise 

camber. For this case the differential equations (A7) become 

^lV) + (VPl") + ( a 3 <P 2 ") + 2,1 ^ 1 ^ 2 ) = P 1 (* 23 ) 

( a 2 <p b") + ( a 3 c Pl" ) + ^4^2 ”) +2p^a 2 cp 2 ^ - 

2(1 - p) + 2^a 2 9 2 'j J = p 2 (A24) 

( a 3 9 0 ") + (a^i") + (a 5 <P 2 ") + Spja^o" + a^" + ^ J - 

4(! - p) |j^8 2 cpi ' J + 2 ^3.^2 ' ) ] + ^ a 1^2 = p 3 (A25) 


with the boundary conditions 


9 0 (°) = <(0) = ^(0) = q> 1 '(0) = cp 2 (0) = q> 2 '(0) = 0 (A2 6) 


/ a ]_9 0 + a 3^2 2 ^ la i c f > 2 j ~ ® 

\ x= 2 

^a 2 cp 0 " + a^" + a^cp 2 " + 2pa 2 qp 2 j = 0 


( a 3 9 0 


+ a]^ + a^qp 2 + Ppa^J = 0 

'x=2 


(A27) 

(A28) 

(A29) 
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[(<W)' + (•2V) 1 + hO' + 2 “( a A)'] l * 0 (*30) 


[kV)’ + ( a 3 ,p l")' + K^")’ + ^ 

2(1 - n) ' + 2 a 2 9 2 ')l = 0 

Jx— Z 

h<)' * (Vi")’ + (Va")' + 

Ml - n)(a 2 <p 1 ' + 2a 3 q) 2 ')J ^ = 0 


(A31) 


(A32) 


If equation (A23) is integrated twice and the boundary condi- 
tions (A27) and (A30) are used, 

„ *2 » a 3 .. 1 r l r l _ , 

<Po =-T- <*>1 ~ a ^2 - 2 ^2 + T / / Pi to2 - (A33) 

a l a l a l J x J x 

/ 

»l « 

Substitution of <Pq into the remaining differential equations and 
boundary conditions results in 


( b l < Pl") + ( b 2 < P2”) ~ 2 ^- ’ + 2(a 2 cp 2 ‘)’] = P 2 " (if J" J Pi ^ 


(A3M 


( b 2 < p 1 ") + (*3^2") - ^ - p) (vO + 2 ( a 3 ( P2 , J' + K 1 ■ p 2 )^ = 


P 3 - 2n 




'a- pZ pZ 


a l J x ^x ?1 


dx2\ " 


(A35) 


( b 1 < P 1 '' + ^2^2") ^ = 0 


(A36) 
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OW) 

1 *1 

+ (b 2 92") - 2(1 - n^a-^' + 2a 2 9 2 ')| = 0 

J x — Z 

(A37) 


(wi' * * 3 * 2 ) - ° 

v x~l 

(A38) 

j(t>2 cp l ) 

+ (b 3 9 2 ") - Ml - |i)(a 2 9 1 ' + 2a 3 9 2 ')l = 0 

-* x=2 

(A39) 


9^0) = 9 1 '(0) = 9 2 (0) = 9 2 '(0) =* 0 

(a4o) 


in which 


2 



If equation (A34) is integrated and the boundary condition (A37) is 

used, 



+ (Va")' • 2(1 





(AM) 


Thus 9^ and 9o are obtained by solving equations (A35) and (AUl) 
with the boundary conditions (A36) , (A38) , (A39) , and (A^O) . Sub- 
sequently, 9 q can be obtained by solving equation (A33) with the 
boundary conditions 9 (0) = 9q'(0) =0. 
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Stresses .- After the approximate deflection of the plate is 
determined from equations (Al8) and (A22) or from equations (A33 ) > 
(A35) , and (A4l) , the extreme-fiber stresses may be calculated from 
the well-known equations of thin-plate theory, which are (see, for 
example, reference 9 ) • 


o 


x 


6d/5 2 w t d 2 w 
t 2 W 2 * By 2 . 


6d /8^w + c&A 
t 2 W 2 Sx 2 / 


T = 6(1 - n)D S 2 w 

xy 2 cbTSy 


The maximum principal stress a at any point in the plate can be 
determined from 
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APPENDIX B 

EXACT SOLUTIONS OF DIFFERENTIAL EQUATIONS FOR 
SOME SPECIFIC DELTA-PLATE PROB LEM S 


The differential equations of appendix A for linear and parabolic 
chordvise deflections are solved exactly for unifor ml y loaded delta 
plates of constant thickness and of diamond chordvise cross section 
vith constant thickness ratio. The equations for deflections obtained 
by the linear theory are presented in terms of the aspect -ratio param- 
eter X for both kinds of delta plates. The equations for deflections 


obtained by' the parabolic theory are presented for — 

1 c 
|i = ^ for the constant-thickness delta plate and for 



and 

h 


\fT 

4 ^ with 
also with 


for the delta plate of diamond chordvise cross section with con- 
stant thickness ratio. 


If the x-axis is passed through the edge perpendicular to the root 
and the substitution x 1 = 1 - j is made, the differential equations 

are clearly of the homogeneous type for which the solutions are of the 
form Xj? where 7 is a constant. For the configurations considered, 
the functions that define the plan form (see fig. l) are then cj(x) = 0 

and C2(x) = cx^, where c is the root chord. In all the equations of 

this appendix the primes denote differentiation with respect to the new 
independent variable x-^. 

Delta Plate of Uniform Thickness under Uniform Load 

Since the stiffness D is a constant for uniform-thickness plates, 
the coefficients in the differential equations (see equation (a 4) ) 
become 



(Bla) 


(Bib) 
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b 2 “ a 4 " 


b 3 " a 5 


■n 


a 2 a 3 Dc 1 * 4 
a x 12 X 1 

(Blc) 

a 3 2 4Dc 5 c 

a x = TT X 1 

(Bid) 

n 

. pi_ x n 

■ n X 1 

(Ble ) 


Solution for linear chordvise deflections .- If the coefficients 
given by equations (Bl) are substituted into equations (A22) and (Al8) 
and the' independent variable is changed to = 1 - the following 

equations for linear chordwise. deflections result ■ 


( x i\") ' " ^ ^i’ ~ -2 x i 3 
V = - | x l ( Pi" + x x 2 


(B2) 

(B3) 


where 

x = ^ 

The boundary conditions to be used with these equations are obtained 
from equations (A13) and (A20) and are 


cp 0 (!) = <P 0 *(1) = 9 X (1) = ^'(l) = 0 . (B4.) 
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The general solution of equation (B2) is 


?! 


= A x 

11 


-1 . -7-1 

+ A 2 x 1 


a! 

4(l - 2X 2 ) Dc 


where 


(B 6) 


7 = + l6X 2 

and A2 and A 2 are arbitrary constants. Since X 2 is inherently 
positive, the boundary condition (B5) requires that A 2 = 0. One 
integration of equation (b 6) and the application of the conditions 
<^(1) = cpi'(l) =0 yields 


<Pl = 


1 pZ 4 

4(1 - 2X 2 ) Dc 


x 7 . 
X 1 



If equation (B3) is solved for <p n with the conditions 
= = 0, the result is u 



1 

1 - 2X 2 






(B7) 


(B8) 


Substitution of equations (B7) and (b 8) into the equation 


w = <p Q + ycp 2 


gives the expression for the deflection w of the plate under the 
assumption of linear chordwise deflections. 

Solution for parabolic chordwise deflections .- If the coefficients 
given by equations (Bl) are substituted into equations (A4l), (A35), 
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and (A33) and the independent variable is again changed to x^ = 1 
the following equations for parabolic chordwise deflections result: 


( X 1 3 V)' + ( X l 4cCp 2")’ " ^(^i' + x i 2cc P 2 ') = " 2 x i 3 ( B 9) 

( x ! V") " + T5( x 1 5c V)" " ^[(x^cp-L*)' + |(x 1 3 ccp 2 ') «j + 


64 4 1 + n 

T . r^-u 





yl k 3 

dT x r 


(BIO) 


Cp 0 " - ~ \ x i c P 1 " - y x^cpg" - 2 hZ 2 cp 2 + ^L. -g- 


(Bll) 


The boundary conditions to be used with these equations are 


9oU) 

= cp 0 '(l) = cPid) = cpi'(l) = <P 2 (l) = 9 2 '(1) 

= 0 

(B12) 


( x i 3 9 1 " + x i 4 ccp 2 ") = 0 


(B13) 


\ ' x 1= 0 




(bV + 15 x i 5c ' p 2") Xi=0 " 0 


(B14) 

( x lV)' 

' + lf( x l 5c<p 2")' - ( x 1 2 cp 1 ' + | x^ccpg') 


= 0 

- 

- 

x =0 



(B15) 
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The homogeneous solutions of the simultaneous equations (B9) 
and (BIO) are of the form 


<Pi' = Ax-, 7 " 1 

-y_l 

<P2 = Bx^ 

Substitution of these expressions into the homogeneous parts of equa- 
tions (B9) and (BIO) leads to the following characteristic equation 
from which 7 may be determined: 


7 6 - 6(l + l6X 2 ) 7^ + |320^4 + + 480X 2 + 9J 

4|i28o + H. \6 + 80^4 + iJ_E^4 + 96x 2 + = o 

and gives the following relationship between A and B: 


7 2 - 


A = -(7 - 1) 


(7 - 2) (7 


(7 + 1) - 1 6 X 2 
- 1 - 16 X 2 


cB 


The particular solutions for uniform loading are given by 

•Pi' = A p x 1 2 
q>2 = Bp*! 2 


(B16) 


where 


A„ = i 


3>i + 1 
1 - u 


- 


2 - 
1 - 


X 2 + 1 


8 - l)x k - (8X2 - l)(4X 2 - l) 


P£ 

Dc 


% " 


1 

¥ 


Vx 2 

8 ^Li(2X 2 - 

l - n v 


( 2 X 2 - l) + 1 + 4x 2 




l)x 4 - ( 8 X 2 - l)(4x 2 - 1 ) DC 2 
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The general solution is the sum of the homogeneous solutions and 
the particular integral 


V = 2_ A n x l 7n_1 + Vl" 

n=l 

V 2 = it. B n x l 7n_1 + Vl 2 


(BIT) 


n=l 


where the values j n are the six roots of the characteristic equa- 
tion (Bl6) and the coefficients A n and B n are the coefficients 
corresponding to each of these roots. After integration cp x becomes 


V - a x l 7n x l^ 

^1 “ 2. An y - + A P ~T~ + A Q 
n=l /r * J 


(B18) 


The general solution for <p Q from equation (Bll) is found to be 


qp 0 


= it. C n x l 7n+1 + C p x l 4 + c q x l + 
n=l 


(B 19 ) 


where, for n = 1, 2, . . .6, 


C n - - 


^n(^n + !) 1 


is— L An + £ 

2 3 




(r„ - i)(r„ - 2 ) ♦ 


B- 




and 


r = - — 
12 


, 2c/, \ 

Ap t ( 1 + — j B P 


The coefficients A-|_ to Ag, Ag, Cg, and C r must be determined by 
the boundary conditions (B12) to (B15). 
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A complete set of coefficients is given in the following table for 
delta plates with Poisson's ratio p. 


equal to 1/3 and with 1 = 1=1 

c 

an<i Deflection curves plotted from these results are shown in 

figures 2 and 3 in which the 45° plate corresponds to 1=1 and the 


a 


60 ° plate corresponds to — = 45 


3 * 


in 

7 1 

n 

A_ 

1 

P z4 

•d Dc 2 
% — T 
pl4 

c -2- 
%z4 

1 = 1 

1 - J2 

3 

1=1 

X.f 

1 = 1 

1 - 42 
3 

1 = 1 

1 _ is 

1 

2 

3 

k 

5 

6 

P 

q 

r 

2.7034 

4.9437 

8.3816 

-2.7034 

-4.9437 

-8.3816 

1.5671 

3.6347 

4.7258 

-1.5671 

-3.6347 

-4.7258 

0.7378 

.02411 

.03827 

0 

0 

0 

-.8000 

-.01557 

0.09632 

.3707 

-.1766 

0 

0 

0 

-.2903 

-.02924 

-0.3133 

-.03039 

-.006293 

0 

0 

0 

.3500 

-0.1022 

-.4313 

.07379 

0 

0 

0 

.4597 

-0.02931 

.003074 

.ooo486 

0 

0 

0 

. o 4 i 67 
-.07152 
.05668 

-0.003223 

.01347 

.002317 

0 

0 

0 

.004032 

-.08354- 

.06692 














Substitution of equations (B17), (Bl8), and (B19) into the equation 


w = cpQ + ytp x + y 2 ^ 

gives the expression for the deflection w of the plate under the 
assumption of parabolic chordwise deflection. 


Delta Plate of Diamond Chordwise Cross Section with 
Constant Thickness Ratio under Uniform Load 


For a delta plate of diamond chordwise cross section with a constant 
thickness ratio the thickness is a function of x and y and is given 
by the following equations: ° 


t = (°SyS^) 

4 - 6 t av( x l - I) S y s cx-l) 
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where t av is the average thickness-. From these expressions for the 
thickness the stiffness can be found and the coefficients in the dif- 
ferential equations become 


a l “ 


a 2 = 


27Dc k 
— X 1 


'n^2 


27Dc 


*i 5 


9Dc3 6 

ao = — — Xn D 




80 




0 X 1‘ 




bp-22si*-7 

2 80 


X l' 


, 26 l 3 Dc^ 8 

3 22400 Xl 


pc n n 

Pn = — *l n 


(B20) 


Solution for linear chordwise deflections .- By use of the coeffi- 
cients given by equations (B20), equations (A22) and (Al8) for linear 
chordwise deflections may be solved for cp-^ and qpQ. The steps in the 

solut'ion are the same in form as those for the uniform- thickness plate 
and the resulting equations are: 


<Pl = 



+ 



(B21) 
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and 

90 " 27 

t lce ^ 1 * 1 ' *1) *£ <B22) 

where 

y = 


Solution for parabolic chordwise deflections .- By use of the coef- 
ficients given by equations (B20), equations (A4l), (A35), and (A33) for 
parabolic chordwise deflections may be solved for 9 ^ q>2' and <?>o* 

The steps in the solution are again the same in form as those for the 
uniform-thickness plate and the resulting general expressions for cp^, 
q> 2 > and cpQ are: 




*1 



\ log e x l 


(B23) 




the exponents 7 are the roots of the characteristic equation 
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OJ zL OJ 

zL l VO 
OJ ^ 



u 

r<e 

OJ 

■H 

it 


CQ 

OJ 


O 


-P 

rH 

0O|CVJ 

a 


<D 

+ 

1 

•H 


O 

1 — 1 1 

Cl 

•H 


N 

ft 

oo| 


ft 


H Zt 

OJ 



+ » 

O 



iH rH 



OJ 


rH O 


t — 

0 

CO CO 

0 

1 _l 

OJ 


0 

^ ^ 

OJ 

CO 

ITAIh. 

cvj 


1 

1 — 1 

ir\U 
cvj r* 

1 

OJ 

It 


Cl 0 
p<ft 

IOJ 

ro\ 0 

pJm 

Olt- 
1 — I ) CVJ 


II 

11 

J 2 - 

ft 


U 

O 

ft 
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C P 2 

A „ 2(1 ^5 qX 2 ' 

1 , . 2 pi 11 
|B P + 51 T 

15 y 

21-Hj 

The coefficients A]_ to 

A 6> A q> 

Cg, and 

C r are again determined by 


the boundary conditions (A26) , (A3 6 ), (A38) , and (A39) in which the coef- 
ficients given by equations (B20) are substituted. 

For Poisson's ratio q equal to 1/3 and X = i = 1, the solution 

of the characteristic equation (B26) leads to two real values and two 
pairs of complex conjugate values for 7. The 'identity 

Xi a± ib H Xl a cos ^ b loggX-j ± ix-^sin^b loggX-j ' 


was therefore used to transform the terms involving the complex conjugate 
values into real form. If - = 1 and q = -k the solution is 

0.004070 Xl 3 - 9Vr - 0.004363 Xl 8 - 075 cos(2.825 log e xJ + 
0.006893 Xl ^’° 7 ^sin(2.825 log eXl j + 0.000294 

-0.003896 Xl 4 ’ 9 ^ 7 + 0.002134 Xi 9 ‘ 075 cos( 2.825 loggXj - 

» 

0.00638l Xl 9,07 ^sin^2.825 logg X;L ) + 0.01794 log^ + 0.001763I 




P z 4 

'pQ “ 


D L 


0.0007715*- 


5.947 


- 0.0000708x- 


10.075 


cos [2 


( 2 

0.001234 Xl 10 ' 075 sin(2.825 loggxj + 0.03331^ 


0.0 4026 


6 _ 


825 loggX-j + • 
log e x l " °* 0l4 - 0 96 Xl 


+ 
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APPENDIX C 


NUMERICAL PROCEDURE FOR SOLVING DIFFERENTIAL EQUATIONS 


In cases where the equations of the present theory cannot he 
solved exactly, a numerical method must he used. In this appendix, 
equations (A19) and equations (A34) and (A35) are set up in difference 
form for numerical solution. Initially the assumption is made that the 
functions involved in the differential equations are continuous and non- 
singular. In this case, first and second derivatives can he expressed 
hy the standard difference forms 

/d 2 y \ _ y n+l ~ 2y n + y n-l 

\dx /„ e 2 



where e is the distance between equally spaced station points. 

In the following development five equally spaced spanwise stations 
are used; however, the extension to a different number of stations can 
he readily made. 

First, consider equation (A19) resulting from the linear theory 

( b i <p i") - 2(1 - ^(Vi')' = *2 - (if Pi<*> - ii 

Because of the nature of the tip boundary conditions for this 
equation, it can he conveniently put in the form 

T' = q 1 (Cl) 

where 

T = ( b i < P 1 ") - 2(1 - n)a 1 <p 1 


f 
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In finding the difference equation equivalent to equation (Cl), the 
quantity ^b-^" j is found in matrix form; from this expression is 

subtracted the matrix equivalent of 2(l - ^a^q^'; the resulting 
expression for . T is multiplied by a differentiating matrix; and the 
product is equated to the right-hand side. 

The quantity the half-stations can be expressed in 

matrix form as follows: 



where the second subscript denotes the station point, the subscript at 
the root station being 0 and at the tip 5* The root boundary conditions 
are now applied; namely. 


q>i(o) = 0 = <p 10 


cp 1 '(o) = 0 


~ ^1-1 

2e 


Thus, after the values of cp iQ = 0 and <p = <p are' substituted, 
equation (C2) becomes 



30 


NACA TN 2621 



One of the tip boundary conditions is 


( b l q ^) X =* ‘ ° " ( blCPl ")5 


Thus, 




The matrix equivalent for the second term of 
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1,7/2 
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The right-hand side of equation (Cl) can now he equated to the 
derivative of equation (C5 ) > thus, 



In order to obtain q.^, the boundary condition 


T = 0 


at x = l must be used. In other words, T goes from sta- 

tion to 0 at station 5. A straight line drawn between these two 

2T9/2 

unnl H Viotro o el nrio iLL 'T'ho value nf* n n c rnnei 


points would have the slope - 
to he this slope; therefore. 


The value of q is considered 
15 




or 
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VO 

o 






■p 03 



<v o 
P <u 
-P rO 



a 

•H 
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In matrix form equation (C7) becomes 
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The foregoing development applies to the case where only linear chordwise deformations are 
allowed. A similar procedure is followed in expressing the differential equations pertaining 
to the parabolic theory in difference form; only the results are shown herein. 



The matrix equivalent to equations (A3*0 and (A35), j 
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,5/2 - ®n,7/2 
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With 9 and cp known^ cp can be obtained by use of equation (A33) 
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It should he noted that, as can be expected, the matrix equa- 
tions (C 6) and (C8) are merely special cases of equations (C9) and (CIO) 
respectively. In addition, the square matrixes in equations (c6) 
and (C9) are symmetric, a result that is consistent with the fact that 
the differential equations under consideration are self-adjoint. 

In the beginning of this appendix the assumption was made that the 
functions involved in the differential equations are continuous and 
nonsingular. The difference solution, however, may be adequate for some 
cases in which this assumption is not strictly correct. For instance, 
the deflections of a plate with a discontinuous stiffness distribution 
could conceivably be not very different from the deflections of a plate 
with a continuous stiffness distribution closely approximating the 
discontinuous distribution except in the neighborhood of the discon- 
tinuity. The results yielded by the difference solution in this case 
would be those associated with the continuous stiffness distribution. 

The number of stations' may have to be increased, however, in order to 
minimize the inaccuracy introduced by the discontinuity or, in other 
cases, by a singularity. The case of the diamond-cross-section, 
constant- thickness -ratio delta plate, discussed in the body of this 
paper, is an example of a treatment of a singularity. In this case, 
although the solution is singular, adequate accuracy is obtained by 
the difference solution if ten equal intervals are used; 
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APPENDIX D 

DEFLECTION AND STRESS EXPERIMENTS ON SOME TRIANGULAR 
CANTILEVER PLATES 


Test specimens .- The specimens tested were: (l) a 1+5° right- 

triangular plate clamped along one leg and (2) a 60° right -triangular 
plate clamped along the longer leg. Each specimen, cut from 
21+S-T1+ aluminum-alloy sheet of 0. 250-inch thickness, had a length per- 
pendicular to the clamped edge of 30 inches. 

Method of testing .- Figure 12, a photograph of the test setup, shows 
the methods of clamping, loading, and measurement of deflections. A 
1,000,000-pound clamping load (held constant during the test) was applied 
to the root area of each specimen and a uniform load of 0.204 psi was 
applied by 2-inch washers giving a tip deflection in each case of approx- 
imately 3/4 inch. 

The deflections were measured by dial gages placed at the points 
indicated in figures 2 and 3* 

Stresses were obtained from the 45° specimen only. On this speci- 
men, 13 resistance-wire rosette strain gages were placed at the points 
indicated in figure 5- The plate was loaded with 2-inch washers in four 
increments of 0.0847 psi per increment and the maximum tip deflection 
was 1.13 inches. Readings of all the strain gages were recorded at each 
increment of loading. 


Analysis and discussion of data . - The deflection w was plotted in 
figures 2 and 3 in terms of the .nondimens ional parameter wD/pZ^, in 

which the elastic constants were taken as E = 10.6 x 10^ psi and p = i. 

It was found that the dial-gage forces reduced the tip deflection of the 
plate by approximately 2 percent) however, since this error is of the 
same order of magnitude as that in the material properties and from 
other sources, no corrections are made in the results presented. 


The readings of each of the 39 individual strain gages were plotted 
against load, and the slope of each of the resulting linear curves was 
taken as the average strain per unit load of the individual gage. The 
principal stresses were then calculated and plotted in figure 5 in terms 
of the nondimensional parameter ot^/pZ^. 
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Figure 2.- Deflections of a delta plate of uniform thickness under 

uniform load. 
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Figure 8.- Numerical and exact . solutions of the differential equations 
for the deflections of the free edges of a 45° delta plate of uniform 
thickness under uniform load. 
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Figure 9*- Numerical and exact solutions of the differential equations 
(obtained by assuming linear chordwise deflections) for the maximum 
principal stresses along the free edges of a 45° delta plate of 
uniform thickness under uniform load. 



re 10.- Numerical and exact solutions of the diff 
obtained by assuming linear chordwise deflections 
long the free edges of a 45° delta plate of diamo: 
ection and constant thickness ratio under uniform 





52 


NACA TN 2621 



Figure 11.- Numerical and exact solutions of the differential equations 
(obtained by assuming linear chordwise deflections) for the maximum 
principal stress along the line y = |(1 - of a b’f delta plate 

of diamond chordwise cross section and constant thickness ratio under 
uniform load. 
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Figure 12.- Deflection test setup of the 45° delta plate under uniform 

load. 


